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Abstract. The solar photospheric carbon abundance has been determined from [Ci], Ci, CH vibration-rotation, CH A-X 
electronic and C2 Swan electronic lines by means of a time-dependent, 3D, hydrodynamical model of the solar atmosphere. 
Departures from LTE have been considered for the C I lines. These turned out to be of increasing importance for stronger lines 
and are crucial to remove a trend in LTE abundances with the strengths of the lines. Very gratifying agreement is found among 
all the atomic and molecular abundance diagnostics in spite of their widely different line formation sensitivities. The mean 
of the solar carbon abundance based on the four primary abundance indicators ([C I], C I, CH vibration-rotation, C2 Swan) is 
log ec = 8.39 ± 0.05, including our best estimate of possible systematic errors. Consistent results also come from the CH 
electronic lines, which we have relegated to a supporting role due to their sensitivity to the line broadening. The new 3D based 
solar C abundance is significantly lower than previously estimated in studies using ID model atmospheres. 
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1. Introduction 

Carbon is a crucial element for the emergence of life in the 
Universe as we know it. Furthermore, it is one of the most 
abundant elements in the cosmos behind only hydrogen, helium 
and oxygen in this respect. While it is clear that carbon is pro- 
duced through the 3a-reaction (Burbidge et al. 1957) in stellar 
interiors, a debate regarding the exact site of its production is 
still ongoing. Recently both low- and intermediate mass stars 
as well as massive stars have been advocated to be the primary 
sources of the element (e.g. Chiappini et al. 2003; Akerman et 
al. 2004). Carbon also plays a key role in the physics of the 
interstellar medium for example through its propensity to form 
dust (e.g. Dopita & Sutherland 2002; Lodders 2003). 

In view of its rather special status in astrophysics, it is 
clearly of particular importance to ascertain the traditional 
abundance anchorage point, namely the solar carbon abun- 
dance. In the standard reference of Anders & Grevesse (1989) 
the solar carbon abundance is quoted as logec = 8.56 ± 
0.04, which was subsequently revised to logec = 8.52 ± 
0.06 by Grevesse & Sauval (1998) based on a re-analysis 
of various atomic and molecular lines using a temperature- 
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modified Holweger-Miiller (1974) ID semi-empirical model 
atmosphere. Holweger (2001) preferred the slightly higher 
value of logec = 8.59 ±0.11 based only on permitted Ci 
lines in the optical and infrared (IR). In sharp contrast to these 
results, Allende Prieto et al. (2002) found the much lower 
value logec = 8.39 ± 0.04 through a detailed study of the 
forbidden [Ci] line at 872.7 nm employing a state-of-the-art 
3D time-dependent hydrodynamical model of the solar pho- 
tosphere (Asplund et al. 2000b). This latter value is in good 
agreement with the evidence from nearby B stars (Sofia & 
Meyer 2001), local interstellar medium (Meyer et al. 1998; 
Andre et al. 2003), the solar corona and the solar wind (Murphy 
et al. 1997; Reames 1999), in particular in view of the recent 
downward revision of the solar oxygen abundance (Allende 
Prieto et al. 2001; Asplund et al. 2004, hereafter Paper IV). 

While the [C I] 872.7 nm line has a well-determined transi- 
tion probability and should be robust against departures from 
LTE, its relative weakness and the possible presence of weak 
blends may cast some doubt on the reliability of the low car- 
bon abundance derived by Allende Prieto et al. (2002). We here 
present a comprehensive study of other available carbon diag- 
nostics, including permitted atomic and various molecular lines 
using the same 3D hydrodynamical solar model atmosphere. 
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Our results corroborate the findings of Allende Prieto et al. and 
unambiguously point to a significant revision by — 0.17dex of 
the solar carbon abundance compared with the often-adopted 
value given in Anders & Grevesse (1989). It is noteworthy that 
only with a 3D analysis are the primary atomic and molec- 
ular abundance criteria in excellent agreement, which further 
strengthens our conclusions. 

2. Analysis 

2. 1. Atomic and molecular data 

[C I] line: Recent calculations of the transition probability of 
the [Cl] 872.7 nm line have apparently converged to A = 
0.640 s" 1 or log gf = -8.136 (Hibbert et al. 1993; Galavis et 
al. 1997). The central wavelength of the line is used as a free pa- 
rameter in the % 2 analysis of the profile fitting described below, 
which results in a rest wavelength of 872.7139 ± 0.0004 nm 
when taking into account the predicted convective blue-shift 
of the line. The lower level excitation potential of the line is 
1.264eV (Bashkin & Stoner 1975). 

Cl lines: The transition probabilities for the selected 16 per- 
mitted C I lines which represent our primary line list were taken 
from Hibbert et al. (1993). The VALD 1 (Piskunov et al. 1995; 
Kupka et al. 1999) and NIST 2 (Wiese et al. 1996) databases 
provided the necessary data for excitation potential, radiative 
broadening and central wavelengths. Collisional broadening 
by H was computed according to the tabulations in Anstee 
& O'Mara (1995), Barklem & O'Mara (1997) and Barklem 
et al. (1998), except for Cl 505.2, 538.0, 658.7, 3085.4 and 
3406.5 nm which fall outside the table boundaries. For those 
lines the classical Unsold (1955) recipe with an enhancement 
factor E = 2.0 was used. The corresponding uncertainty is 
however minor: with E = 1.0 the derived abundances for the 
three optical lines would be about 0.01 dex higher while for the 
two IR lines the corresponding error is about 0.03 dex. 

We initially considered additional Cl lines, which how- 
ever were all rejected in the final analysis for various reasons, 
such as suffering from blends (e.g. 477.5, 1 186.3 nm), too large 
sensitivity to velocity and/or pressure broadening (e.g. 965.8, 
1068.5, 1070.7, 1189.5, 1744.8, 3372.8nm) or uncertain con- 
tinuum placement (e.g. 1 184.8 nm). We stress that it is better 
to retain only the best quality abundance indicators than using 
a larger sample of less reliable lines. Nevertheless, our line list 
is sufficiently large to secure a final mean abundance that is 
insensitive to possible problems with a few individual lines. 
CH vibration-rotation and electronic lines: The analysis is 
based on the best 102 vibration-rotation lines of the (1,0), (2,1) 
and (3,2) bands in the infrared (3.3 — 3.8 ^m). We also include, 
with a lower weight, nine apparently clean weak lines from the 
(0,0) and (1,1) bands of CH A-X around 430 nm (6 and 3 lines, 
respectively). The dissociation energy of CH is well-known: 
A)(CH) = 3.465 eV (Huber & Herzberg 1979). The parti- 
tion functions and statistical weights for CH are taken from 
Sauval & Tatum (1984) for consistency with the adopted disso- 
ciation energy and equilibrium constants. The lower level ex- 
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citation potentials for the CH lines were adopted from Melen 
et al. (1989) and the transition probabilities from Follmeg et 
al. (1987). The latter are the same as used by Grevesse et al. 
(1991). The band oscillator strengths for the CH A-X lines are 
/oo = (5.16 ± 0.11) • 10~ 3 and / n = (4.46 ± 0.15) • 10" 3 
(Larsson & Siegbahn 1986). As the quantum mechanical ap- 
proach for pressure broadening outlined by O'Mara and col- 
laborators (e.g. Anstee & O'Mara 1995) is not applicable to 
molecular transitions, we have resorted to the classical Unsold 
(1955) recipe with an enhancement factor of 2 for all molec- 
ular lines. The particular choice of enhancement factor has no 
impact on the derived abundances however for these lines. 
C 2 electronic lines: We selected the 17 least blended weak C 2 
lines from the (0,0) d 3 U g - a 3 Tl u Swan band. Our adopted 
dissociation energy of C 2 comes from the accurate experimen- 
tal measurement by Urdahl et al. (1991): D (C 2 ) = 6.297eV. 
This laboratory determination is in very good agreement with 
the theoretical calculations of Pradham et al. (1994) which sug- 
gests £> (C 2 ) = 6.27 eV. Our value is 0.09 eV higher than that 
recommended by Huber & Herzberg ( 1 979) and 0. 1 9 eV higher 
than the value employed by Lambert (1978) in his classical so- 
lar abundance analysis of CNO. We remind the reader that the 
derived carbon abundance from C 2 lines depends on the disso- 
ciation energy as Alogec ~ —ADo/2. The C 2 partition func- 
tions and equilibrium constants have been adopted from Sauval 
& Tatum (1984). The C 2 Swan band oscillator strengths are 
well-determined from measurements of the radiative lifetime 
of the d 3 Tl g state (see Grevesse et al. 1991). As for the CH 
lines, we use an enhancement factor of 2 to the Unsold recipe 
for pressure broadening, which however has no influence on 
the final results. We do not consider the lines of the C 2 Phillips 
system to be sufficiently reliable for our purposes and hence do 
not include any such lines in the analysis. 
CO vibration-rotation lines: In principle, CO lines can be 
utilised to derive the carbon abundance whenever the oxygen 
abundance has first been determined accurately from alterna- 
tive diagnostics. We will not consider CO in this analysis, how- 
ever, due to their large model atmosphere sensitivity through 
the high dissociation energy which makes them less reliable 
than other diagnostics. Indeed, no available ID hydrostatic 
model atmosphere, be it semi-empirical like the Holweger- 
Miiller (1974) or VAL-3C (Vernazza et al. 1976) models or 
theoretical like the MARCS (Gustafsson et al. 1975 and sub- 
sequent updates) or Kurucz (1993) models, have been able to 
explain the strong CO lines observed in the Sun (not count- 
ing of course semi-empirical models specifically constructed 
to reproduce the CO lines, such as the models described by 
Grevesse & Sauval 1991, 1994, Grevesse et al. 1995 and Avrett 
1995). This has been taken as evidence for a very cool tem- 
perature structure or evidence for a temperature bifurcation 
of the gas in the upper photosphere/lower chromosphere (e.g. 
Ayres 2002 and references therein). Uitenbroek (2000a,b) has 
shown that the situation has improved with the advent of ID 
time-dependent hydrodynamical atmosphere models which ac- 
count for important non-LTE effects (Carlsson & Stein 1992, 
1995, 1997) and 3D hydrodynamical models like those de- 
scribed herein. Indeed, with the same 3D model as used here 
the long-standing COnundrum may finally be about to be re- 
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solved (Asensio Ramos et al. 2003). We intend to return to the 
CO lines in a future investigation in this series. 

2.2. Observational data 

Throughout we use disk-center (/i = 1.0 for the optical 
lines and /i = 0.935 for molecular lines in the IR) intensity 
observations. In the optical wavelength region the Brault & 
Neckel (1987, see also Neckel 1999) solar atlas has been em- 
ployed. This high-quality atlas was recorded with the Fourier 
Transform Spectrograph at Kitt Peak Observatory at a resolv- 
ing power of about 500,000 and a signal-to-noise exceeding 
1000 at all wavelength regions of interest here. Minor adjust- 
ments to the tabulated continuum level immediately surround- 
ing the relevant spectral lines were made. In some cases the 
Kitt Peak FTS atlas was supplemented with the Jungfraujoch 
disk-center solar atlas (Delbouille et al. 1973) when lines were 
affected by telluric lines (e.g. C I 71 1.1 and 960.3 nm). For the 
four lines between 1254 and 1259 nm we used the infrared 
solar atlas of Delbouille et al. (1981) also recorded at Kitt 
Peak. For the analysis of the last four Cl lines in Tableland 
the CH vibration-rotation lines the Spacelab-3 ATMOS 3 solar 
disk-center intensity IR observations has been used (Farmer & 
Norton 1989, see also Farmer 1994). The resolving power of 
the ATMOS atlas is about 200,000. The S/N varies with wave- 
length but is typically at least 400 in the regions of interest 
here. 

A major advantage with the new generation of 3D hydrody- 
namical model atmospheres employed here is the in general ex- 
cellent agreement between predicted and observed line profiles 
without resorting to the usual micro- and macroturbulence nec- 
essary in ID spectrum synthesis (e.g. Asplund et al. 2000a,b,c). 
This facilitates the use of line profile fitting for individual lines 
for the carbon abundance determination, a method that we ap- 
ply here to the analysis of the atomic lines. Since the line profile 
fitting procedure is more ambiguous in the ID case due to the 
poorer agreement between predicted and observed profiles, we 
here adopt the theoretical line intensities from the 3D profile fit- 
ting as "observed" equivalent widths for the ID spectrum syn- 
thesis (we note that these values are in excellent agreement with 
those directly measured on the observed spectra). This enables 
a direct comparison of the effects stemming from the choice 
of ID or 3D model atmospheres. Due to the sheer number of 
molecular lines employed in the present study, we rely for prac- 
tical reasons on the measured equivalent widths for both the ID 
and 3D abundance analyses of the molecular transitions. The 
equivalent widths of the molecular lines have been obtained by 
fitting Gaussians to the observed profiles. Various tests have 
ensured that no significant error has been introduced by the use 
of fitting Gaussians rather than for example a Lorentz profile 
for these particular lines. 

2.3. LTE spectral line formation 

We follow the same procedure as in previous articles in the 
present series (Asplund et al. 2000b,c; Asplund 2000; Asplund 
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et al. 2004; Asplund 2004; hereafter Papers I, II, III, IV, V) of 
studies of spectral line formation in the solar granulation. We 
employ a 3D, time-dependent, hydrodynamical simulation of 
the solar surface convection as a realistic model of the solar at- 
mosphere. The reader is referred to Stein & Nordlund (1998) 
and Paper I for further details of the construction of the 3D solar 
model atmosphere. For comparison purposes, we have also per- 
formed identical calculations employing two widely-used ID 
hydrostatic models of the solar atmosphere: the semi-empirical 
Holweger-Miiller (1974) model and a theoretical, LTE, line- 
blanketed MARCS model (Asplund et al. 1997). 

Equipped with the 3D hydrodynamical solar atmosphere, 
3D spectral line formation calculations have been performed 
for [C I], C I, CH, and C2 lines under the simplifying assump- 
tions of LTE. In addition, instantaneous chemical equilibrium 
(ICE) has been assumed valid for the molecule formation. 
This is a potentially serious issue since in a stellar atmosphere 
the time-scale for molecule formation could be longer than 
the corresponding dynamical time-scale and hence the actual 
molecular number densities may be out of equilibrium (e.g. 
Uitenbroek 2000a,b). To our knowledge, the only published 
study devoted to follow the time-dependent chemical evolution 
using a reaction network is that of Asensio Ramos et al. (2003) 
who investigated CO in the Sun. They concluded that ICE is 
an acceptable approximation below heights of about 700 km. 
Since our CH and C2 lines are formed in much deeper lay- 
ers, this assumption is unlikely to significantly affect our de- 
rived C abundances (see also Asensio Ramos & Trujillo Bueno 
2003). However, we stress that only detailed chemical evolu- 
tion calculations can confirm whether this conclusion is cor- 
rect. While such computations are beyond the scope of the 
present study, we intend to investigate this in a forthcoming 
publication. Finally we note for completeness that in principle 
it is necessary to use the non-LTE populations of H and C when 
estimating the number densities of CH and C2 even under the 
assumption of ICE. Fortunately, both elements are predomi- 
nantly populated by their ground states for which LTE is an ex- 
cellent assumption in the relevant parts of the solar atmosphere. 
Therefore the use of the Saha distribution when computing the 
atomic partial pressures needed for estimating the molecular 
densities is fully justified. 

The spectrum synthesis employs realistic equation-of-state 
and continuous opacities (Gustafsson et al. 1975 and subse- 
quent updates; Mihalas et al. 1988). The temporal coverage of 
the part of the solar simulation which forms the basis of the 
line formation computations is 50 min of solar-time with snap- 
shots every 30 s. The comparison with observations relates to 
disk-center intensity profiles (/1 = 1.0 and /1 = 0.935 de- 
pending on the line). As the calculations self-consistently ac- 
count for the Doppler shifts arising from the convective mo- 
tions, no micro- or macroturbulence enter the 3D spectral syn- 
thesis. In the absence of such convective line broadening, all ID 
calculations have been performed adopting a microturbulence 
fturb = l.Okms -1 unless noted otherwise (e.g. Holweger & 
Miiller 1974: Blackwell et al. 1995). 
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2.4. Non-LTE spectral line formation 

While non-LTE effects are expected to be modest for the other 
carbon abundance diagnostics, the same can not be assumed 
for the high-excitation permitted C I lines employed here. In or- 
der to investigate this further we have performed detailed non- 
LTE calculations based on the Holweger-Miiller and MARCS 
ID model atmospheres. Full 3D non-LTE line formation com- 
putations have recently been carried out for Li (Asplund et al. 
2003) and O (Paper IV) for a couple of snapshots from the same 
solar simulation as employed here. As described below, the re- 
quirement for a very extensive carbon model atom to obtain 
realistic results currently prevents similar calculations for this 
element. Until such a 3D non-LTE study is performed we will 
have to rely on the results from the corresponding ID case, 
which, however, is expected to yield sufficiently accurate re- 
sults for these high-excitation lines formed in deep atmospheric 
layers. In this connection, we note the very similar non-LTE 
abundance corrections obtained in 3D and ID for Ol lines of 
similar character as the C I lines (Paper IV; Allende Prieto et 
al. 2004). We intend to return to the issue of 3D non-LTE line 
formation for C I lines in a later study. 

Our carbon model atom consists of 217 levels in total, of 
which 207 belong to C I, nine to C II and one represents C III. 
The atom is complete up to principal quantum number n = 9 
for the singlet and triplet systems with no fine structure of the 
terms accounted for. The highest C I level corresponds to an ex- 
citation potential Xcxc = 11.13eV, which is only 0.13eV be- 
low the ionization limit. The energy levels are taken from the 
experimental data of Bashkin & Stoner (1975). The levels are 
coupled radiatively through 458 bound-bound and 216 photo- 
ionization transitions with the data coming from the Opacity 
Project database TOPbase (Cunto et al. 1993). Excitation and 
ionization due to collisions with electrons are accounted for 
through the impact approximation for radiatively allowed tran- 
sitions and the van Regemorter (1962) formalism for forbid- 
den lines adopting an oscillator strength / = 0.01. No in- 
elastic hydrogen collisions are included as the available labora- 
tory measurements and quantum mechanical calculations sug- 
gest that the normally employed recipe (Drawin 1968, see also 
Steenbock & Holweger 1984) over-estimates the collisional 
cross-sections by about three orders of magnitude (Fleck et al. 
1991; Belyaev et al. 1999; Barklem et al. 2003). We do note, 
however, that there is some astrophysical evidence that this is 
not the case for O (e.g. Allende Prieto et al. 2004). In the ab- 
sence of suggestions to the contrary for C, we prefer not to em- 
ploy the Drawin formula given the uncertainties involved. The 
ID non-LTE calculations have been performed using the statis- 
tical equilibrium code MULTI (Carlsson 1986), version 2.2. The 
non-LTE abundance corrections have been interpolated from 
the non-LTE and LTE equivalent widths for four different abun- 
dances (logec = 8.30,8.40,8.50,8.60). The corresponding 
LTE results necessary for estimating the non-LTE abundance 
corrections were also obtained with MULTI using a model atom 
with all collisional cross-sections artificially set to extremely 
large values to ensure full thermalization of all levels and tran- 
sitions. 
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Fig. 1. The departure coefficients (3[ — rii/np™ for the Cl 
levels as a function of continuum optical depth in the MARCS 
model atmosphere; the corresponding case for the Holweger- 
Miiller model is almost identical. The four lowest levels (dotted 
lines) show no departures from LTE, while the following two 
levels (dashed lines) experience a slight over-population in the 
line-forming region. The majority of levels, however, have a 
small under-population (solid lines), as does the ground state 
of C II (solid line connected with stars). 



Table [2 lists the resulting non-LTE abundance corrections 
for the Holweger-Miiller and MARCS ID model atmospheres. 
Clearly, the non-LTE effects are not particularly dependent on 
the employed model atmosphere, at least not within the frame- 
work of ID models. While the four lowest C I levels (2p 2 3 P, 
2p 2 1 D, 2p 2 1 S, 2p 3 5 S°) are all perfectly described by the 
Boltzmann distribution due to close collisional coupling, the 
next two levels (3s 3 P°, 3s 1 P°) are slightly over-populated rel- 
ative to the LTE expectation, as seen in Fig.^ The vast major- 
ity of levels, however, show a very similar behaviour with small 
but significant underpopulations in the typical line-forming re- 
gions (—1 <> logT5oo ^ according to the line depres- 
sion contribution functions introduced by Magain 1986): /3; = 
n-JnY 1 ^ < 1. The similarity of the levels with x cxc > 8.5eV 
is mainly due to efficient collisional interlocking. In all cases, 
the non-LTE abundance corrections for the here employed C I 
lines are negative (i.e. the lines become stronger in non-LTE 
than in LTE), since the line source function exceeds the Planck 
function S\/B v k, /3 up //3i ow < 1 in the line-forming regions. 
In addition, many of the Cl lines originate from the 3s 3 P° 
or 3s 1 P° levels which are slightly over-populated (J3i ow > 1) 
which increases the line opacity and hence strengthens the lines 
further. The abundance corrections show a clear dependence 
on equivalent width with the largest corrections obtained for 
the strongest lines. Our non-LTE abundance corrections would 
have been significantly smaller (< 0.05 dex) had we included 
H collisions through the classical Drawin (1968) formula. Our 
non-LTE results are qualitatively the same as those obtained 
by Sturenbock & Holweger (1990) but show in general slightly 
larger non-LTE effects, mainly due to our larger atom, more 
realistic photo-ionization cross-sections and no inclusion of ef- 
ficient H collisions. 



M. Asplund et al.: Solar line formation: VI. The photospheric C abundance 



5 



Table 1. The computed non-LTE abundance corrections of 
Cl lines for the Holweger-Muller (1974) model and MARCS 
(Asplund et al. 1997) ID model atmospheres. The departures 
from LTE have been estimated using MULTI (Carlsson 1986) 
and a 217-level model atom. 



line 
[nm] 



Xexc 

[eV] 



logs/ 



A (logec) 



HM 



MARCS 



505.2167 

538.0337 

658.7610 

711.1469 

711.3179 

960.3036 

1075.3976 

1177.7546 

1254.9493 

1256.2124 

1256.9042 

1258.1585 

2102.3151 

2290.6565 

3085.4621 

3406.5790 



3s 1 ?? - 

3P 1 ?! ■ 

3p 3 D! 
3p 3 D 3 
3s 3 Pg - 
3s 3 P§ - 
3p 3 D 2 
3p 3 Po - 
3p 3 Pi - 
3p 3 Pi - 
3p 3 Pi - 
3p 1 S - 
3p 1 S 
Sd^ 
4p x Pi 



4p 1 D 2 

^d'P? 
- 4d 3 F 2 ' 
- 4d 3 F2 
-3p 3 Si 

- 3p 3 Di 

- 3d 3 F§ 

- 3d 3 P1 

- 3d 3 Pg 

- 3d 3 P? 

- 3d 3 Pa 

- Sd^? 
-4s 1 P? 
-4p 1 D 2 
-4d x D§ 



7.685 
7.685 
8.537 
8.640 
8.647 
7.480 
7.488 
8.643 
8.847 



9.172 
9.172 
9.736 
9.989 



-1.304 
-1.615 
-1.021 
-1.074 
-0.762 
-0.895 
-1.598 
-0.490 
-0.545 
-0.504 
-0.586 
-0.509 
-0.437 
-0.182 
0.086 
0.454 



-0.05 
-0.04 
-0.03 
-0.04 
-0.05 
-0.15 
-0.11 
-0.11 
-0.09 
-0.09 
-0.09 
-0.09 
-0.08 
-0.10 
-0.05 
-0.11 



-0.04 
-0.03 
-0.03 
-0.04 
-0.05 
-0.15 
-0.10 
-0.09 
-0.08 
-0.08 
-0.08 
-0.08 
-0.06 
-0.08 
-0.05 
-0.10 



3. The solar photospheric C abundance 

3. 1. The forbidden [Ci] 872. 7nm line 

The forbidden [C I] 872.7 nm line has been studied in detail by 
Allende Prieto et al. (2002) using the same 3D hydrodynami- 
cal solar model atmosphere as employed here. They determined 
the carbon abundance from a x 2 -analysis of the solar flux atlas 
of Kurucz et al. (1984) with three free parameters besides the 
carbon abundance: the wavelength of the [C I] line (see discus- 
sion in Sect. 12. 1> . continuum placement and log gf ■ eg; for the 
Si I 872.80 nm line of which the [C I] line is located in its blue 
wing. Including fitting errors and estimates of possible system- 
atic errors, they arrived at logec = 8.39 ± 0.04. 

For consistency with the analysis of the permitted lines, we 
have re-analysed the forbidden line using the disk-center in- 
tensity profile instead of the flux profile. The derived carbon 
abundance, however, is not affected by this choice: logec = 
8.39 ± 0.04. We note that Lambert & Swings (1967) flagged 
a potential blending Fe I line, which, if contributing to the fea- 
ture, would imply a downward revision of our derived carbon 
abundances. Using other Fe I lines belonging to the same multi- 
plet, Allende Prieto et al. (2002) estimated the equivalent width 
of the line to be < 0.03 pm (1 pm = 10 mA), which corresponds 
to a downward change of < 0.02 dex for the derived carbon 
abundance. 

The best fit carbon abundance results in a theoretical disk- 
center intensity equivalent width of 0.53 pm. Employing this 
value, the ID Holweger-Muller model atmosphere implies a 
carbon abundance of log ec = 8.45 ± 0.04 (adopting the same 
error estimate as for the 3D analysis), while the MARCS model 
suggests logec = 8.40±0.04. Clearly, this particular transition 
is relatively insensitive to the choice of model atmosphere. 



Table 2. The derived solar carbon abundance as indicated by 
the forbidden [C I] 872.7 nm transition and permitted C I lines. 
The results include the non-LTE abundance corrections pre- 
sented in Tabled In the absence of detailed 3D non-LTE cal- 
culations for Cl, the ID non-LTE corrections for the MARCS 
model atmosphere has here been used also for the 3D case. 



line 






3D a 




HM C 


MARCS C 


[nm] 


[eV] 






[pm] 






[Cl]: 














872.71 


1.264 


-8.136 


8.39 


0.53 


8.45 


8.40 


Cl: 














505.2167 


7.685 


-1.304 


8.34 


4.07 


8.40 


8.34 


538.0337 


7.685 


-1.615 


8.36 


2.52 


8.42 
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8.33 
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8.34 
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8.35 


5.53 


8.36 


8.34 
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9.989 


0.454 


8.34 


7.09 


8.34 


8.34 



The abundances derived using the 3D model atmosphere have 
been obtained from profile fitting of the observed lines. 
The predicted disk-center intensity equivalent widths with the 3D 
model atmosphere using the best fit abundances shown in the 
fourth column. 

The derived abundances with the Holweger-Muller (1974) and the 
MARCS (Asplund et al. 1997) ID hydrostatic model atmospheres 
in order to reproduce the equivalent widths presented in the fifth 
column. 



3.2. Permitted Ci lines 

The carbon abundance determination using permitted atomic 
lines is based on the 16 C I lines listed in Table [2] The indi- 
vidual abundances have been estimated using profile fitting for 
disk-center intensity; we note for completeness that very simi- 
lar abundances would have been derived had flux profiles been 
used instead. The profile fitting has been performed for the spa- 
tially and temporally averaged 3D LTE line calculations from 
the 100 snapshots used. As clear from Fig. |2 the agreement 
between predicted and observed profiles is very satisfactory. 
In contrast, the theoretical line profiles in ID require additional 
broadening from macroturbulence to get the correct line widths 
but of course the observed line shifts and asymmetries can not 
be accounted for. 

Due to the large carbon model atom required to capture 
the essence of the non-LTE effects, 3D non-LTE line forma- 
tion computations for the C I lines have not been possible. As 
a substitute, we apply the ID non-LTE abundance corrections 
estimated with the MARCS model atmosphere given in Table 
^to the 3D LTE abundance estimates to arrive at our final re- 
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Fig. 2. The predicted 3D LTE profiles of the C I 71 1.3, 1075.4, 1 177.7 and 2290.6 nm lines (diamonds) compared with the solar 
intensity atlas (solid lines, Brault & Neckel 1987 for the first three lines and Farmer & Norton 1989 for the last line). The excellent 
agreement is achieved without any micro- or macroturbulence due to the Doppler shifts arising from the convective motions. 



suits presented in Table It is worthwhile pointing out that 
the similarities between the ID non-LTE results for the MARCS 
and Holweger-Miiller model atmospheres, in spite of their quite 
different temperature structures (see Fig. 2 in Paper IV), give 
reasons to believe that the 3D non-LTE case would indeed be 
similar as well had it been available. 

The mean C 1-based solar carbon non-LTE abundance is 
logec = 8.36 ± 0.03 (standard deviation). The scatter is en- 
couragingly small. There is no trend in the derived abundances 
with either wavelength or excitation potential although as seen 
in Fig.[3]there is a small correlation with the strengths of the 
lines (about 0.04 dex from the weakest to the strongest lines, 
i.e. W\ = 1 — 12 pm). One possibility is that the line broad- 
ening arising from the Doppler shifts due to the convective 
motions are underestimated, akin to using too small a micro- 
turbulence in ID. Our previous 3D LTE analysis of Fe I lines 
have revealed a similar trend (Asplund et al. 2000c). However 
it should be pointed out that that correlation could also be the 
result of neglecting departures from LTE (Shchukina & Trujillo 
Bueno 2001), and there is no corresponding trend for Fe II or 
Si I lines (Asplund 2000), which would be expected if the con- 
vective line broadening is underestimated. From Fig. [3] it is 
clear that the trend with equivalent width is largely driven by 
the two weakest lines, 658.7 and 711.1 nm. Although the for- 



mer line is located in the far wing of the Ha line and there are 
some telluric H2O lines around the latter, we have not been able 
to identify any obvious explanation for these two abundances 
being underestimated. 

In our opinion, the most likely explanation for this slight 
correlation is underestimated non-LTE corrections in the 3D 
case. We note that we have here adopted the non-LTE calcu- 
lations from the ID MARCS model in the lack of suitable full 
3D non-LTE calculations. It is quite conceivable that the pres- 
ence of temperature inhomogeneities could produce slightly 
more pronounced non-LTE effects, in particular for the stronger 
lines. Fortunately, this residual trend is very small and will 
not affect the mean abundance from the C I lines significantly. 
Another possibility, while less likely, is that the ^/-values 
for the Cl 658.7 and 711.1nm are over-estimated by about 
0.05 dex. 

We emphasize that the trend with equivalent width would 
have been more pronounced, as seen in Fig. [3] had departures 
from LTE not been considered (Fig. [5} and the line-to-line 
scatter would have been significantly higher (±0.05 instead of 
±0.03). This is a clear indication that these highly excited C I 
lines are indeed not formed in LTE. 

Using the equivalent widths estimated from the 3D profile 
fitting and the available non-LTE corrections, the ID carbon 
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abundances become log ec = 8.35 ± 0.03 with the MARCS and 
logec = 8.39 ± 0.03 with the Holweger-Miiller model atmo- 
spheres. As for the 3D case, without application of the com- 
puted non-LTE corrections there would be a pronounced trend 
with equivalent widths (Fig. [3}- However, in both ID cases the 
non-LTE abundances exhibit a trend in the opposite sense with 
respect to the LTE case. 

3.3. CH vibration -rotation lines 

High-quality, high resolution solar IR spectra like the Spacelab- 
3 ATMOS atlas open new opportunities for element abundance 
determinations compared with the traditional UV-optical re- 
gion employed in most studies. One of the major advantages 
with the IR is the relatively clean spectrum with few blending 
lines. In addition, many noteworthy molecules have their vi- 
bration and rotation lines accessible in this region, which can 
provide a multitude of unperturbed lines to base an abundance 
analysis on. Grevesse et al. (1991) reported the first solar car- 
bon abundance determination using CH vibration-rotation lines 
in the IR and argued that these lines constitute one of the most 
reliable sources for determinations of the solar carbon abun- 
dance. Here we agree with this conclusion, but note that this 
is only true when the lines are analysed with a realistic model 
of the solar atmosphere. Due to the temperature sensitivity of 
molecule formation in general (e.g. Asplund & Garcia Perez 
2001), the strengths of molecular transitions and consequently 
the derived abundances depend crucially on the photospheric 
temperature structure. In addition, the presence of atmospheric 
inhomogeneities induced by for example the convective mo- 
tions can also strongly alter the equivalent widths of the lines. 
As these concerns are properly addressed by our use of a real- 
istic 3D hydrodynamical solar model atmosphere, we consider 
the CH vibration-rotation lines as primary abundance indica- 
tors. 

In this study we rely on 102 CH vibration-rotation lines 
of the (1,0), (2,1) and (3,2) bands around 3.3-3.8 ^m, for 
which the equivalent widths were measured in the Spacelab- 
3 ATMOS IR intensity (p = 0.935, i.e. nearly at solar disk- 
center) atlas. With our 3D model and assuming LTE for the line 
formation and molecular concentration, these CH lines imply a 
carbon abundance of logec = 8.38 ± 0.04. The derived abun- 
dances for each line together with the employed line data are 
listed in Table [3] As seen in Fig.|4j there are no trends in de- 
rived carbon abundance with wavelength or equivalent width, 
although there is a weak correlation with lower excitation po- 
tential (0.03 ± O.Oldex/eV). It should be noted however that 
the total span in excitation potential is only 1 .4 e V. We do not 
consider this weak trend as flagging a serious problem, in par- 
ticular given the encouragingly small scatter. These lines are 
also formed in a rather narrow region of the atmosphere (mean 
optical depth of line formation —1.25 logrsoo —1.15 in 
the Holweger-Miiller model atmosphere; line formation depth 
is a ill-defined concept in a 3D model atmosphere). We have 
not been able to identify any particular problem with the CH 
3527.3 nm line, which gives a significantly lower C abundance 
(log ec = 8.25) than the other lines. 



The corresponding result using the Holweger-Miiller model 
atmosphere gives a substantially higher result: logec = 8.53 ± 
0.04. This is a direct consequence of the higher temperatures in 
this ID model. In this case the abundance trend with excitation 
potential has disappeared and instead been replaced with a mi- 
nor correlation with equivalent width. This opposite behaviour 
compared with the 3D case reflects the correlation between 
equivalent width and excitation potential of the employed CH 
lines. With the MARCS model atmosphere, the CH lines indi- 
cate a carbon abundance of log ec = 8.42 ± 0.04 without any 
apparent trends with transition properties. 

3.4. C 2 electronic lines 

There are numerous weak lines from the C2 Swan band in the 
solar spectrum. From these we selected a subsample of 17 (0,0) 
lines between 495 and 516nm, which are all apparently undis- 
turbed by neighboring lines. The lines employed here have 
equivalent widths between 0.3 and 1.4 pm, which make the de- 
rived abundances insensitive to the velocity broadening. Since 
the previous uncertainties surrounding the C2 dissociation en- 
ergy appear to have dissolved (see Sect. 12. 1> . we include these 
lines among our primary abundance indicators. 

With the 3D hydrodynamical solar model atmosphere the 
measured equivalent widths of these 17 C2 lines yield a so- 
lar carbon abundance of logec = 8.44 ± 0.03. There are 
no significant abundance trends with line properties accord- 
ing to the results presented in Table [4] The C2 based abun- 
dance is in excellent agreement with the values derived from 
the other preferred diagnostics. The corresponding results for 
the ID models are loge c = 8.53 ± 0.03 (Holweger-Miiller) 
and log e c = 8.46 ± 0.03 (marcs). 

3.5. CH electronic lines 

While there is a swath of CH electronic lines in the optical solar 
spectrum the vast majority are badly blended with other lines or 
are too strong to yield accurate results. We have identified only 
nine apparently unblended lines from the (0,0) and (1,1) bands 
of CH A-X between 421 and 436 nm. Their equivalent widths 
are between 3.5 and 8 pm, making them partly saturated and 
thus sensitive to the atmospheric velocity field. That together 
with the rather crowded spectral region the lines are located in, 
relegates the CH electronic lines to carry only a secondary role 
in our solar C abundance analysis. 

The measured equivalent widths of these nine CH lines to- 
gether with the 3D hydrodynamical model atmosphere yield a 
C abundance of logec = 8.45 ± 0.04. No significant trends 
with wavelength, excitation potential or equivalent width are 
present (Table |5}, although the spans in these parameters are 
quite small. The derived 3D abundance is slightly higher than 
those from the primary abundance indicators (Cl, [Cl], CH 
vibration and C2 electronic lines). The corresponding results 
with the Holweger-Miiller and MARCS model atmospheres are 
loge c = 8.59 ± 0.04 and loge c = 8.44 ± 0.04, respec- 
tively. As for the CH vibration lines, the higher temperatures in 
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Table 3. The derived solar carbon abundance from CH vibration-rotation lines. The individual abundances have been derived 
using measured disk-center (p = 0.935) intensity equivalent widths for both the 3D model as well as the Holweger-Miiller and 
MARCS ID model atmospheres. 
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the Holweger-Miiller model automatically lead to higher abun- 3.6. Summary 
dances. 

Of the five carbon abundance indicators considered here, the 
greatest weight is given to the forbidden [Cl], permitted Cl, 
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Table 4. The derived solar carbon abundance as indicated by 
lines from the C2 (0,0) Swan band. As for the other carbon 
abundance diagnostics the abundances have been derived for 
the 3D model and the Holweger-Miiller and MARCS model at- 
mospheres. The measured equivalent widths are for disk center 
intensity (p = 1.0). 



line Xcxc l°g#/ W\ 3D HM MARCS 

[nm] [eV] [pm] 



495.1400 


1.070 


0.631 


0.50 


8.46 


8.55 


8.47 


499.2300 


0.820 


0.572 


0.65 


8.43 


8.51 


8.44 


503.3700 


0.570 


0.181 


0.46 


8.42 


8.51 


8.43 


503.7700 


0.580 


0.486 


1.15 


8.48 


8.57 


8.50 


505.2600 


0.500 


0.450 


1.00 


8.42 


8.51 


8.44 


507.3600 


0.370 


0.073 


0.56 


8.41 


8.50 


8.43 


508.6200 


0.340 


0.345 


1.35 


8.47 


8.56 


8.49 


510.3700 


0.260 


0.267 


1.25 


8.45 


8.54 


8.47 


510.9100 


0.230 


0.237 


1.25 


8.45 


8.54 


8.47 


510.9300 


0.220 


-0.088 


0.54 


8.41 


8.50 


8.43 


513.2500 


0.140 


-0.243 


0.49 


8.42 


8.52 


8.44 


513.6600 


0.110 


-0.343 


0.53 


8.47 


8.57 


8.50 


514.0400 


0.100 


-0.405 


0.37 


8.42 


8.52 


8.44 


514.3300 


0.100 


-0.414 


0.46 


8.47 


8.57 


8.49 


514.4900 


0.100 


-0.450 


0.34 


8.42 


8.52 


8.44 


515.0500 


0.340 


0.332 


1.03 


8.41 


8.50 


8.43 


516.3400 


0.170 


0.113 


0.93 


8.40 


8.50 


8.43 



Table 5. The derived solar carbon abundance based on CH A- 
X electronic lines in the optical using three different model at- 
mospheres. The listed equivalent widths are measured in the 
disk-center (/i = 1.0) intensity solar atlas. 



line 


Xexc 


log gf 


W X 


3D 


HM 


MARCS 


[nm] 


[eV] 




[pm] 








421.8723 


0.413 


-1.008 


7.88 


8.41 


8.55 


8.40 


424.8945 


0.192 


-1.423 


6.52 


8.42 


8.56 


8.41 


425.3003 


0.523 


-1.523 


3.70 


8.44 


8.57 


8.43 


425.3209 


0.523 


-1.486 


4.00 


8.45 


8.58 


8.44 


425.5252 


0.157 


-1.453 


6.80 


8.45 


8.59 


8.44 


426.3976 


0.460 


-1.593 


3.90 


8.48 


8.61 


8.47 


427.4186 


0.074 


-1.558 


6.59 


8.45 


8.59 


8.44 


435.6375 


0.157 


-1.830 


4.74 


8.54 


8.68 


8.53 


435.6600 


0.157 


-1.775 


4.20 


8.40 


8.54 


8.40 



CH vibration-rotation and C2 electronic lines, with the CH 
electronic lines only assigned a supporting role in view of 
their broadening sensitivity and location in a relatively crowded 
spectral region. The four primary abundance indicators imply 
highly concordant results when analysed with the 3D hydro- 
dynamical model atmosphere, as summarised in Table [6] This 
is particularly noteworthy given the very different tempera- 
ture sensitivity and distinct formation depths of these lines. In 
sharp contrast, the Holweger-Miiller model atmosphere gives 
much more disparate results, with the molecular lines suggest- 
ing much higher abundances than the atomic transitions as a 
consequence of the different temperature structures and lack of 
photospheric inhomogeneities. The theoretical MARCS model, 
however, performs nearly as well as the 3D model in this re- 



Table 6. The derived solar carbon abundance as indicated by 
forbidden [C I], permitted C I, CH vibration-rotation, C2 elec- 
tronic and CH electronic lines. The results for the C I lines in- 
clude non-LTE abundance corrections. The quoted uncertain- 
ties only reflect the line-to-line scatter for the different types of 
C diagnostics, except for [C I] where the error estimate includes 
uncertainties in the profile fitting. 



lines 




logec 






3D 


HM 


MARCS 


primary indicators: 








[CI] 


8.39 ±0.04 


8.45 ± 0.04 


8.40 ±0.04 


Ci 


8.36 ±0.03 


8.39 ± 0.03 


8.35 ±0.03 


CH vib-rot 


8.38 ±0.04 


8.53 ± 0.04 


8.42 ± 0.04 


C2 electronic 


8.44 ±0.03 


8.53 ± 0.03 


8.46 ±0.03 


secondary indicator: 








CH electronic 


8.45 ±0.04 


8.59 ± 0.04 


8.44 ±0.04 



spect; a similar conclusion would likely have held had we used 
for example a Kurucz theoretical ID model atmosphere. 

The errors quoted in Table |6] for the different types of 
transitions reflect the line-to-line scatter rather than smaller 
standard deviation of the mean, except for the forbidden [C I] 
872.7 nm which also includes estimates for fitting errors and 
systematic errors. The mean of the primary indicators becomes 
logec = 8.39 ± 0.03. The dominant source of error though 
is likely to be of systematic nature. The excellent agreement 
between the different types of lines gives us some confidence 
that unknown systematic errors can not be very significant and 
estimate them to be on the order of ±0.05 for the mean abun- 
dance. We therefore arrive at our best estimate of the solar car- 
bon abundance as: 

logec = 8.39 ±0.05. 



4. Comparison with previous studies 

The solar carbon abundance derived in Sect. [3] is significantly 
lower than most previous estimates. In his extensive and thor- 
ough survey of the solar C, N and O abundances, Lambert 
(1978) found logec = 8.69 ± 0.10 based on a combination 
of [Ci], CH A-X and C 2 lines using the Holweger-Miiller 
semi-empirical model atmosphere. He preferred to exclude the 
permitted C I lines from the final average due to uncertainties 
stemming from the available transition probabilities and pos- 
sible departures from LTE, which in the time since then have 
at least partly been addressed. Likewise, the CH B-X, CH C-X 
and CO vibration lines, neither of which have been considered 
here, were relegated to a supporting role in Lambert's study for 
a variety of reasons. He did not have access to the CH vibra- 
tion lines which we rank as a primary abundance indicator. In 
contrast, we consider the CH A-X lines as inferior to the other 
indicators. The main differences with Lambert's value for [C I] 
compared with our own re-analysis using the Holweger-Miiller 
model originate in the adopted gj-values (— 0.08dex) and pre- 
ferred disk-center equivalent width (— O.lOdex) but we have 
not been able to trace the remaining — 0.06dex discrepancy. 
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The use of a 3D model atmosphere instead of the Holweger- 
Miiller model further decreases the derived carbon abundance 
by 0.06 dex. 

Lambert's carefully derived C abundance was subsequently 
duly adopted by Grevesse (1984) in his review of the solar 
abundances. Soon thereafter, Sauval & Grevesse (1985) identi- 
fied CH vibration lines in the solar IR spectrum, which should 
enable a reliable abundance determination, as also demon- 
strated herein. In their widely used survey of the solar and 
meteoritic abundances, Anders & Grevesse (1989) preferred 
logec = 8.56 ± 0.04 based primarily on a preliminary anal- 
ysis of the CH vibration lines using Holweger-Miiller model 
atmosphere. Once the C analysis was complete, this value had 
been slightly revised to logec = 8.60 ± 0.05 based on 104 
CH vibration lines (Grevesse et al. 1991). In addition, Grevesse 
et al. analysed C2 (Swan and Phillips bands), CH A-X and C I 
lines and estimated a mean solar carbon abundance of log ec = 
8.60 ± 0.05. The main reason for the high molecular based 
abundances compared with ours clearly stems from the use of 
the Holweger-Miiller model instead of a 3D hydrodynamical 
model atmosphere, which takes into account temperature inho- 
mogeneities aand has a cooler mean temperature stratification 
(Table |6j. The difference with their C I abundances is primar- 
ily a combination of them neglecting non-LTE effects, choice 
of ^/-values and their adopted equivalent widths, with a minor 
part probably coming from the employed pressure broadening 
data. 

More recently, Grevesse & Sauval (1999) performed a re- 
analysis of the C I, [C I], CH and C2 lines. Based on the exis- 
tence of significant trends in the obtained Fe abundances with 
excitation potential with the standard Holweger-Miiller model 
atmosphere, Grevesse et al. (1999) derived a new temperature 
structure which was slightly cooler in the higher layers. This 
modified Holweger-Miiller model atmosphere forms the basis 
for the slightly lower solar C abundance estimated by Grevesse 
& Sauval (1998): loge c = 8.52 ± 0.06. As before, the ne- 
glect of photospheric inhomogeneities causes the molecular 
lines still to yield a too high abundance. 

As described above, Allende Prieto et al. (2002) performed 
a detailed study of the [C I] using the same 3D hydrodynam- 
ical solar model atmosphere as employed in the present work 
and found a much lower value than other studies: logec = 
8.39 ± 0.04. This low value is corroborated here by our more 
extensive calculations also involving permitted C I and molec- 
ular lines. In contrast, Holweger (2001) concluded that the 
solar carbon abundance is logec = 8.59 ± 0.11 based on 
Cl and [Cl] lines when applying non-LTE (—0.05 dex) and 
granulation corrections (+0.02 dex) to the ID LTE results of 
Sturenburg & Holweger (1990) and Biemont et al. (1993). We 
note that these granulation corrections have the opposite sign 
to those estimated here (Table |2j as a result of Holweger's 
different definition of these (Paper IV). The main differences 
to our Holweger-Miiller-based results lie in the selection of 
lines (we restricted our list largely to lines with W\ ^ 10 pm 
while Holweger's sample contains a large number of lines with 
10 ^ W\ Si 20 pm), choice of ^/-values, adopted equiva- 
lent widths and the inclusion or not of inelastic H collisions in 



the non-LTE calculations. Holweger did not consider molecu- 
lar lines. 

In summary, we are confident that our comprehensive 3D 
analysis superseeds previous works devoted to the solar photo- 
spheric carbon abundance. This conclusion is supported among 
other things by the excellent agreement between the different 
abundance indicators illustrated in Table[6] We emphasize that 
the transition from a ID hydrostatic to a 3D hydrodynamical 
model atmosphere is only part of the explanation for the large 
downward revision of the solar carbon abundance compared 
with for example Lambert (1978) and Grevesse et al. (1991). 
Improved atomic and molecular data, higher quality observa- 
tions, accounting for non-LTE effects and identification of ex- 
isting blends also played important roles. 



5. Conclusions 

We have presented a comprehensive study of the available 
atomic and molecular lines to derive a new solar carbon abun- 
dance which is significantly lower than most previous analy- 
ses: log e c = 8.39 ± 0.05 (C/H = 245 ± 30 • lO" 6 ). The quoted 
uncertainty includes our best estimate of the remaining system- 
atic errors. Our primary abundance indicators are the forbidden 
[Cl] 872.7 nm, high-excitation permitted Cl, CH vibration- 
rotation and C2 Swan lines with a supporting role played by CH 
A-X electronic lines. A special effort has been made to select 
the most accurate and reliable atomic and molecular data such 
as transition probabilities and dissociation energies. For the C I 
lines detailed non-LTE calculations in a ID context have been 
performed, which revealed significant departures from LTE in 
the line formation. Finally, a novel feature of the present study 
is the use of a realistic 3D hydrodynamical solar model atmo- 
sphere, which furthermore has proved to be crucial in order to 
obtain consistent abundances between the various diagnostics. 
In sharp contrast, the ID Holweger-Miiller model atmosphere 
yields much higher abundances for the molecular lines than for 
the atomic lines, partly due to the higher temperatures in the 
line-forming regions and partly due to the neglect of tempera- 
ture inhomogeneities. The excellent concordance between the 
various transitions with widely different temperature and pres- 
sure sensitivities is a very strong argument in favour of both the 
3D model atmosphere as such and the new low carbon abun- 
dance. 

Further support for the new low solar carbon abundance 
advocated herein comes from a comparison of the chemical 
composition in related environments. The revised solar pho- 
tospheric carbon abundance is now in good agreement with 
those measured in nearby B stars with solar iron abundances: 
logec = 8.28 ±0.17 or C/H = 190 ± 90 • 10~ 6 (Sofia & Meyer 
2001). Furthermore, the solar carbon abundance is now in good 
agreement with estimates of the total (gas plus dust) local inter- 
stellar medium carbon abundance for realistic gas-to-dust ratios 
(Andre et al. 2003; Estiban et al. 2004). Finally, coupled with 
the previously re-determined solar photospheric oxygen abun- 
dance (Paper IV), the photospheric C/O ratio is 0.54 ± 0.09, 
which is in good agreement with the measurements in solar 
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flares (0.42 ± 0.09, Fludra et al. 1999; 0.54 ± 0.04, Murphy et 
al. 1997) and solar wind particles (0.47±0.01, Reames 1999) 4 . 

The new solar carbon abundance corresponds to a change 
of — 0.28dex relative to that advocated in Lambert (1978) and 
— 0.17dex compared with the preferred value in the standard 
source of Anders & Grevesse (1989). This change is primarily 
the outcome of improved models of the solar photosphere and 
the line formation processes, as well as more accurate atomic 
and molecular line data. The possibility of similar substantial 
systematic errors in standard ID abundance analyses still today 
is certainly worth keeping in mind when interpreting the results 
of other abundance studies in terms of stellar nucleosynthesis 
and galactic evolution. 
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Fig. 3. The derived solar carbon abundance from permitted C I 
lines using the 3D hydrodynamical model atmosphere as a 
function of equivalent width (Upper panel). The filled circles 
denote the case when the ID MARCS-based non-LTE abun- 
dance corrections in Table^have been applied, while the open 
circles correspond to the LTE results. Also shown are the corre- 
sponding cases with the Holweger-Miiller (Middle panel) and 
the MARCS (Lower panel) ID model atmospheres, which dis- 
play similar trends with equivalent widths. As explained in the 
text, this may signal underestimated non-LTE effects. 
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Fig. 4. The derived solar carbon abundance (filled circles) from 
CH vibration-rotation lines using the 3D hydrodynamical time- 
dependent simulation of the solar atmosphere (Asplund et al. 
2000b) as a function of wavelength, lower level excitation po- 
tential and equivalent width (in pm). The solid lines denote 
least-square-fits giving equal weights to all lines. 



